Entanglement entropy in one-dimensional disordered interacting system: The role of 

localization 
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The properties of tiie entanglement entropy (EE) in one-dimensional disordered interacting sys- 
tems are studied. Anderson localization leaves a clear signature on the average EE, as it saturates 
on length scale exceeding the localization length. This is verified by numerically calculating the EE 
for an ensemble of disordered realizations using density matrix renormalization group (DMRG). A 
heuristic expression describing the dependence of the EE on the localization length, which takes into 
account finite size effects, is proposed. This is used to extract the localization length as function of 
the interaction strength. The localization length dependence on the interaction fits nicely with the 
expectations. 
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Introduction. — Although more than half a century has 
passed since the seminal work of Anderson on the absence 
of diffusion in certain random lattices [l|, the localiza- 
tion transition, with an emphasis on understanding its 
behavior in strongly correlated systems, remains one of 
the central themes of condensed matter physics. Unlike 
many other quantum phase transitions in which a gap 
between the ground state and the excited state energies 
appear, the localization transition is more subtle. The 
canonical manifestation of the localized phase is the ex- 
ponential dependence of the conductance, G, on the lin- 
ear dimension of the system, L. The localization length, 
^, (corresponding to the correlation length in the local- 
ized regime) is defined through the exponential decrease 
in the conductance G{L) ^ exp(— i/^) y . 

For non-interacting systems the localization length can 
be visualized easily as the the length-scale on which the 
single-electrons wave-function envelope decays. No such 
transparent interpretation exists for the many-particle 
wave-function. Indeed, the numerical extraction of the 
localization length for interacting systems is difficult. Di- 
rect calculation of the conductance is computationally 
taxing since it demands knowledge of the excited states. 
Ground state properties are easier to calculate, but there 
one runs into other problems. For example, using the 
sensitivity to boundary conditions [3| measured, e.g., by 
the persistent current, depends both on the localization 
length as well as on the inverse compressibility of the sys- 
tem [J]. Calculating the two requires both a calculation 
of the sensitivity of the ground state to a magnetic flux in 
periodic boundary conditions as well as the dependence 
of the number of electrons in the system on the chemical 
potential. Alternatively, one may study the response of 
the system to the introduction of a perturbation [5| . For 
example, the decay of the Friedel oscillations resulting 
from the introduction of an additional impurity to the 
system [6|. Nevertheless, this is strictly speaking justi- 
fied only for weak disorder, and one must calculate the 
ground state both in the presence and absence of the 
impurity. 



In this letter we will utilize the entanglement entropy 
[3] (EE, sometimes referred to as the von Neumann 
or Shanon entropy) in order to extract the localization 
length of a ID interacting disordered system. EE for a 
many-body system in a pure state is related to parti- 
tioning it into two regions: A and B. The entanglement 
between A and B is measured by the EE Sa/b in the fol- 
lowing way: Using the Schmidt decomposition the many- 
body pure state of the entire system, |^), is expressed 
as the sum of two orthonormal basis sets of regions A 
({|0A,i)}) and B {{\4>b,])}), in the following way 



I*) ^T.^a,\(|)A,^)®\<i)B,^), 



(1) 



with S^a^ — 1. 



The Schmidt decomposition is related to the eigenbasis 
of the reduced density operators Pa/b = Tr5/^|^)(4'|, 
by 



Pa/b = ^i^i \4>A/B,i){4'A/i 



(2) 



The EE measures the entanglement between the two re- 
gions by the von Neuman entropy of the reduced density 
matrix: 



Sa/b = -Sia^ln(a^), 



(3) 



equivalent to the Shannon entropy of the squared 
Schmidt coefficients af , clearly Sa — Sb- 

The behavior of the EE has recently been connected 
to the presence of quantum phase transitions (QPTs) in 
condensed mater systems [7H10|. Specifically, the EE of 
a region of length La of a ID metallic system is expected 
to grow logarithmically with La [llMl3j , while the EE in 
the insulating regime should not depend on the regions 
size. Since for length scales shorter than the localization 
length the system behaves essentially as a metal, we ex- 
pect that the EE of a disordered ID system will show a 
logarithmic dependence for La ^ £,, while it will saturate 
for La ^ C- Tracing the crossover will be used to deter- 
mine the localization length ^, even for an many-particle 
interacting ground state. 



Density matrix renormalization group (DMRG) [IJ, 
|l5| is the natural choice for a numerical method to cal- 
culate the ground state properties of disordered interact- 
ing ID system. Moreover, since the Schmidt coefficients 
are calculated in the DMRG in order to decide which 
states should be truncated, calculating S has no addi- 
tional computational overhead. Since DMRG is much 
more accurate for finite systems with open boundaries 
than for periodic boundary conditions (another drawback 
for methods relaying on persistent current calculations), 
we prefer to consider a finite system of size L in which 
the EE of a region of length La starting at the edge is 
calculated. Thus, the EE will depend both on La as well 
as on L. The explicit function S{La, L) as been derived 
by several authors Ill4l3lll6l|. and for open boundaries: 
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S{La-, L) = - In ( sin 



f^LA\ 






(4) 



where c is a non-universal constant which depends also 
on L and the boundary entropy 
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This dependence 
was tested numerically for ID spin chains [18|, and found 
to fit pretty well. 

Model. — We consider a spinless ID electrons wire of 
size L with repulsive nearest neighbor (NN) interactions 
and on-site disorder. The system's Hamiltonian is thus 
given by the Anderson model 

L L-l 

L-l 



^YA""! ~ 9)(4+l'^J + l - x), 



where e^ is the random on-site energy, taken from a uni- 
form distribution in the range [—W/2, W/2], U is the NN 
interaction strength ([/ > 0), and i = 1 is the hopping 
matrix element between NN. cj is the creation operator 
of a spinless electron at site j in the wire, and a positive 
background is included in the interaction term. For the 
non-interacting case the dependence of the localization 
length on the disorder is known ^{W, U = 0) ~ 105/W^ 

[19.]. 

For the interacting case, using renormalization group 
[20J the following dependence of the localization length 
on interaction strength and disorder is suggested 



^(W^, [/) = (e(W^, [/ = 0))l/(3-29(C/))^ 



(6) 



where q(U) = 7r/[2 cos^^(— [//2)] is the Luttinger param- 
eter [2l|. For non-interacting electrons g{U = 0) = 1. 
Since for repulsive interactions g < 1 decreases as a func- 
tion of the interaction strength, one finds that the lo- 
calization length always decreases as a function of the 
interaction strength 



22| 



Clean systems. — We begin by calculating the EE for 
clean systems of different sizes as function of the length 
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FIG. 1: (Color online) The EE of a clean system as function 
of La, for different system sizes L. The symbols correspond 
to the DMRG results, the dotted curves to Eq. (U) and the 
continuous curves to Eq. (O, with a value of ("o indicated 
in the legend. Bottom inset: The EE of a clean system for 
different values of interaction strength U and L — 700. The 
symbols correspond to the DMRG results, the dotted curves 
to Eq. Q and the continuous curves to Eq. (O, with (o = 
255.6. Top inset: i^o as function of L. The curve corresponds 
to Co(i) = 0.357L. 



of region A (La). During the DMRG calculations 192 
states are retained. The accuracy may be assessed by 
comparing S{La,L) to S{L — La,L). The results are 
presented in Fig. [TJ The EE for all system length begin 
with a logarithmic dependence on La for La <C L and 
tapers out for La ~ L/2. Fitting the numerical results to 
Eq. (U) , where the only fitting parameter is the constant 
c, works quite well, except for some deviations at small 
values of La- Such deviations were also seen for spin 
chains with open boundary conditions (see Ref. llSl ). 

Adding interactions {U 7^ 0) to a clean system does not 
change the conductance and the system remains metallic. 
Thus, one expects that the EE will behave in a similar 
way in the presence or absence of interactions. Indeed, 
this can be seen in the bottom inset of Fig. [TJ where 
the EE as function of La is calculated for different in- 
teraction strength {U — 0,0.15,0.3,0.6 corresponding to 
g = 1,0.954,0.913,0.836), for L = 700. Except for a 
small change in the constant c, the EE shows the same 
functional behavior for the different values of interaction. 

Disorder. — The study of the infiuence of disorder on 
the EE of ID systems has concentrated on the infiuence 
of a single impurity in the middle of the sam ple 
|29|. For the non-interacting Anderson model [29 
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the 

system remains metallic even in the presence of an impu- 
rity and the EE continues to be described by a logarithm, 
though with a different slope. For the interacting case, 
(i.e. for a Luttinger liquid) the impurity effectively cuts 
the system into two, driving the system into an insulat- 
ing phase [30| , resulting in a saturation of the EE [23| . 



The multiple impurity situation was only discussed for 
the non-interacting case [3l|, |32| , where the EE of a sin- 
gle cite (i.e., La = I) was shown to correspond to the 
well known measure of the inverse participation ratio. 

When disorder is introduced, the main changes antic- 
ipated in the behavior of S{La,L) are the saturation of 
of the EE for L > La ^ £, and a metallic behavior 
(i.e., logarithmic, however with a non- universal slope) 
for i,^ ^ La- The following considerations are taken 
into account while suggesting a functional description of 
the EE: For small values of x = L^/min{^, L/2} one ex- 
pects S{La, Ci L) k, S{x) k, \n{x) while for large values of 
X, S{x) ~ Constant. Thus, it is natural to assume that 
S{x) takes the form S{x) « ln(/(a;)), where f{x) « x for 
a; ^ 1 and f{x) w Constant for a; ^ 1. The most natural 
choice would be f{x) = tanh(a;), nevertheless it turns out 
that the crossover between the two regions {f{x) w x and 
f{x) « Constant) is quite sharp and tanh(a;) does not de- 
scribe it well. The straight forward way to sharpen the 
transition would be to replace the exponential governing 
tanh(a;) by a stronger decaying function, e.g. a Gaussian, 
hence f{x) = erf(a;) (the error function) which provides 
a better description. Thus, we propose to fit the EE to 
the following heuristic expression: 



5(L^,f,L)-s(01n erf 



La 



(7) 



where s(^) < 1/6, is a constant which depends on the 
disorder, ({L, ^) depends both on L and ^ and c is a non- 
universal constant. It is important to note that for the 
disordered case one discusses a disorder averaged value 
of the EE. 

For the clean limit (^ — ;> cx)), we may infer the behavior 
of Co{L) = C(i, oo) from fitting Eq. ([7]) to the data we 
obtained for different system length presented in Fig. [T] 
The fit is represented by the continuous curves in Fig. [T] 
As expected s(^ — >■ oo) « 1/6, while the results for C,o{L) 
are depicted in the legend. One can see that the heuristic 
expression for the EE (Eq. ([7])) provides a very good 
description of the numerical data with an even somewhat 
better fit than Eq. (U) for small values of La- In the 
top inset of Fig. [1] Co as function of L is presented. A 
linear relation Co = 0.357L emerges. Eq. (O fits also the 
behavior of the EE for the clean interacting case (bottom 
inset Fig. [1]) where only c depends on the interaction 
strength. 

In the limit of C *C i one would expect that the sys- 
tem size L will have only a rudimentary influence on 
C(i,^), i.e., C,{L » 5) « ^. We investigate this asser- 
tion by studying the behavior of S{La, C, L) for W = 0.7 
and U = 0.6 for L > 300. Using Eq. ([5]), one expects 
that the localization length ^ — 57, is much smaller than 
L. Prior to examining the behavior of the average value 
S{La,^t L), it is important to examine the behavior of 
the distribution of the EE for different realizations, in or- 
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FIG. 2: (Color online) The median of the EE {S{La,^,L)) 
as function of La for system length L = 300, 500, 700. The 
median was calculated for ensembles of 600, 400, 400 corre- 
spondingly. In all cases, the disorder W — 0.7 (^ = 214) 
and the interaction strength U = 0.6 {g = 0.836). Error bar 
represent range between the 40th and 60th percentiles. The 
curves represent a fit to Eq. ((Tjl with values of C indicated 
in the legend. Left inset: The probability distribution of the 
EE for L = 700 and La = 20, 330. The symbols represent the 
numerical data, while the curves correspond to a fit to the 
normal distribution. Right inset: The same numerical data 
as in the left inset, fitted to a log-normal distribution. 



der to verify that the concept of an average value of the 
EE is meaningful, especially since we are in the localized 
regime. The distribution of S{La,£,,L) for L — 700 de- 
rived from an ensemble of 400 realizations and different 
values of La is presented in the insets of Fig. [2] In the 
left inset the distribution is fitted to a normal distribu- 
tion while in the right inset to a log-normal distribution. 
While for La = 20 the distribution follows a normal dis- 
tribution for La = 330 a fit to the log-normal distribution 
works better. Thus, it seems, that while for La ^ £, the 
EE distribution is normal, for La 3> ^ it is log-normal, 
somewhat similar to the behavior of the distribution of 
the conductance [33|. Thus, a better representation of 
the typical EE is the median, which is less sensitive to 
the tails of the distribution. The errors are estimated by 
the 40 (and 60) percentile of the distribution. 

The median EE for samples of length L = 300, 500, 700 
over ensembles of 600,400,400 realizations correspond- 
ingly are presented in Fig. [21 It is clear that the EE 
does not change much for the different length of the sam- 
ples. Fitting Eq. ^ results in C = 63.6.69.6,63 for 
L = 300, 500, 700. This is in line with our expectations 
for C » ^ for L > C- 

For the case where both the finite-size of the sample 
and the localization length are of comparable length, one 
expects that both will influence C(i,C)- Since the limits 
are clear, i.e., C{L > ^) w C and C{L < C) = Co = 
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FIG. 3: (Color online) As in Fig. [2l for weak interactions 
([/ = 0.15,5 = 0.954)). The curves represent a fit to Eq. ^ 
with values of ^ indicated in the legend. Inset: Comparison 
of the interpolation formula proposed in Eq. ^ represented 
by the curve, and the numerical results for ^ (circles). 



0.357L, one may interpolate between these two limits by: 
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''0 
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(8) 



In Fig. [3]we examine the behavior of the EE for a weakly 
interacting {U — 0.15, g = 0.954) case. According to Eq. 
([S]), ^ = 136 is close to the values of Co expected for the 
system sizes L = 300, 500, 700 examined. A fit of the 
numerical data to Eq. ^ resuhs in C = 79.9, 92.3, 120.7 
for L = 300, 500, 700. As can be seen in the inset, these 
values of C fit reasonably with the extrapolation of C ob- 
tained from Eq. (j5]). 

Now we examine the role played by interactions on the 
behavior of the localization length. In Fig. U] we show 
the EE of the longest samples (L = 700) for which it 
was feasible to obtain an adequate number of realiza- 
tions. As expected the saturation of the EE begins for 
smaller values of La as the interaction increases. A more 
quantitative comparison of the behavior of the localiza- 
tion length may be achieved by fitting the median EE to 
Eq. ^ and extracting C,. Finite size corrections (espe- 
cially important when C, is comparable to (^Oj i-e. for weak 
interactions) can be taken into account by the extrapo- 
lation suggested in Eq. (|S|) . The localization length thus 
obtained is compared to the expectations according to 
Eq. (O in the inset. We see that the localization length 
obtained from the EE fits quite well the expectations for 
the dependence of ^ on the interaction strength. 

Conclusions. — We have studied the properties of the 
EE in one-dimensional disordered interacting systems. 
The EE saturates once the length of the segment exceeds 
the localization length. We propose a heuristic expres- 
sion which describes the dependence of the EE on the 
segment length, system size and localization length and 




FIG. 4: (Color online) The median of the EE as function 
of La for system length L = 700 and W = 0.7 (^ = 214) 
for different values of interaction strength U = Q, 0.15, 0.3, 0.6 
corresponding to g = 1, 0.954, 0.913, 0.836. The curves repre- 
sent a fit to Eq. ((Tj with values of (" indicated in the legend. 
Inset: The values of C, (depicted by circles) and ^ (extracted 
from Eq. ^ using C,o = 255.6, indicated by squares) as func- 
tion of the Luttinger parameter g. The curve corresponds to 
the values of ^ expected according to Eq. ©. 



compare it to numerical results. By fitting the expression 
to the numerical data it is possible to extract the local- 
ization length of an interacting system while calculating 
only its ground state many-body wavefunction. 
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